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Xsokine^lc Sampling 

The present invention relates to isokinetic 
sampling. It is particularly, but not exclusively, 
5 related to methods and apparatuses for caxrrying out 
isokinetic sampling, and in particular isokinetic 
sampling in oilfield applications. 

In oilfield applications, as in many other fields, 
10 it is important to be able to "analyse the coii5>osition and 
properties of a multi-phase fluid stream, such as a gas- 
condensate stream or a wet-gas stream. In particular it 
is desired to know the quantity of the various phases of 
the flow, for exair^le by knowing the gas^oil ratio (GOR) 
15 or the condensate-gas ratio (CGR) . It is also desired to 
• know the properties of the various phases, for exan^le 
their pressure, volume and ten5>erature (PVT) 
relationships . 

For example, the prediction of fluid phase behavior 
20 and reservoir simulation models based on equations-of- 
state (EOS) typically requires high quality PVT data. 

PVT data are also often needed for flow assurance in 
wells and transport lines . 

Furthermore, PVT data often have a significant 
25 impact on processing facility designs and specifications, 
and therefore on the profitability of gas-oil field. 

However, to accurately determine the PVT properties 
and coitposition of reservoir fluid, representative fluid 
sartroles are required. Common analysis sampling procedures 
30 suffer from deficiencies in either or both of their 

accuracy and their ability to cope with high flow rates. 
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For example, a high-rate gas condensate well 
producing 5-100 MMscf (Millions standard cubic feet) /day 
with a gas-oil ratio of 3k-100k scf/bbl {CGR from 0 to 
200 bbl/MMscf) (bbl = barrel of liquid; 1 = 6.29 bbl) , 
5 can often overwhelm a test-separator, causing liquid 
carry-over- in the separator gas outlet line, thereby 
providing poor measurement of the GOR and non- 
representative PVT samples (if the carry-over is not 
measured and corrected for) . This normally also results. 

10 in poor recombination ratios. Gas condensate wells are 
particularly problematic for sampling as large volumes 
of gas are associated with only small volumes of liquid, 
and the phase behavior will be highly sensitive to the 
quantity and composition of the liquid phase. 

15 Wellhead sampling is therefore considered to be the 

only practical method of obtaining reliable test data for 
such applications. Isokinetic sampling at a wellhead of 
a multiphase fluid is desirable since, if achieved, it 
means that the sample of the fluid is at the same 

20 pressure, temperature and velocity as the main flow 

stream, and therefore will have identical properties to 
that main flow stream. 

A known wellhead sampling system is shown 
25 schematically in Figure 1, which is taken from Dybdahl 
and Hijermstad ^'A systematic approach to sampling during 
well testing". Society of Petroleum Engineers (SPE) paper 
69.427, SPE Latin America & Caribbean Petroleum 
Engineering Conference, Buenos, Argentina, 25-28 March 
-30 2001 (see also the paper by Henk et al. '"Testing of gas 
condensate reservoirs - sampling, test design and 
analysis", SPE paper 68668, SPE Asia. Pacific Oil & Gas 
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Conference, Jakarta, Indonesia, 17-19 April 2001). Such 
a system- is, for example, produced by Petrotech under 
the±r IsoSplit(TM) wellhead sampling and separator 
sarnpling services . 

The wellhead sampling system^ of Figure 1 consists of 
two parts . 

The first part is a wellhead mixing/sampling 
manifold 1010, which is positioned upstream of the choke 
man±fold 103 0, and has a f ocusing-action static mixing 
device 1015 (see GB 2301297 & US 5S94080; a similar 
mixing device is shown in Figure 2a) for collecting 
liquid at the wall of the pipe and distributing it into 
the main flow stream. A saitpling probe 1020 is radially 
inserted downstream of the mixer 1015 for withdrawing a 
samfDle of the main flow stream. 

The second part is an on-site laboratory 1050 (also 
called a ^^Mini-lab" ) , which processes the sampled stream 
to measure the condensate-gas ratio (CGR) . PVT samples 
for compositional characterization are taken. 

As described in SPE paper 68668, the liquid 
distribution across the flow pipe of Figure 1 will vary, 
and either a mixing device or a traversing probe has to • 
be used in order to have a split stream sait^le with the 
corr-ect gas condensate ratio. As described in SPE paper 
6942 7, the probe traversing enables sampling at multiple 
poiiits across the pipe and the subsequent correction for 
the droplet and velocity distribution. It is stated in 
SPE 68668 that the liquid collected from the wall by the 
focu sing-action mixer 1015 is distributed homogeneously 
into the flowing gas stream. However, our own gas-liquid 
ratio (GLR) profile measiirement shown in Figure 2b, for a 
mixear design (shown in Figure 2a) similar to the 
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focussing-action mixer 1015^ indicates a liquid-rich jet 
at the . central region of the sairpling cross section- The 
GLR measured at such points will be a grross under- 
estimate of the pipe average value. Figure 2b shows a - 
7 0% error in the GLR at the pipe centre, with respect to 
tlae pipe reference GLR. 

SPE paper 68658 also states that sample withdrawal 
tlarough the sampling probe 20 from the main stream is 
performed at an isokinetic sampling rate, i.e. at the 
same . Linear velocity as the wellhead fluid stream, which 
is necessary to sample the gas and liquid phase in the 
correct ratio. This isokinetic sampling rate is said to 
be calculated from the probe size and the gas flow rate 
measured at the test separator. 

This approach has at least the following 
dx sadvantages : 

1) - A test separator is needed for wellhead 

sampling operations, which is bulky and 
expensive (Figure 1) . 

2) The test separator can be overwhelmed by 
high-rate gaig wells, leading to liquid 
carry-over in the separator gas outlet line. 
This can result in an inaccurate main-stream 
gas rate measurement used to set the 
isokinetic sampling rate. 

.3) The main-flow' gas rate can fluctuate, for 

example due to cheuiges in the liquid 
loajding. This \insteadiness in the main flow 
rate can :cause the separator gas rate 
measurement to be \mrelis±>le (see SPE paper 
76766, by B.C. Theuveny and P. Mehdizadeh 
^*Multiphase flowmeter application for well 
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and fiscal allocation", Alaska, 20-22 May 
2002), leading to xinreliable isokinetic 
ScLinpling. 

Isokinetic sait^ling rate control based on 



measured main-stream gas flow rates can lead 
to errors when liquid loading becomes 
relatively high. 
WOOO/49370 describes a method for measuring the 
liqu±d and gas flow rates of a multi-phase fluid stream 
inclTJiding withdrawing a portion of the overall flow rate 
undenr ""isokinetic conditions". The sampled portion forms 
approximately 5% to 15% of the total flow rate. 

The methods and apparatus disclosed in WOOO/49370 do 
not ±ndependently measure the main flow stream, but 
instead derive measurements of the main flow stream from . 
the effect of the sampling operation on the flow stream. 
Furtliermore , although the methods described are claimed 
to be isokinetic", this term is applied to all sainpling 
in which the ratio of the sampled flow rate to the 
calcTjilated total flow rate is within 20% of the ratio of 
the oross sectional area of the sampling probe to the 
cross sectional area of the flow stream. Adjustment of 
the sampling conditions is only made if the ratio falls 
outside this range. 

The method proposed in wbOO/49370 for measuring the 
total flow rate (Qi) of the fluid stream is given in 
equation (5) of that application, with Qx being calculated 



is tile measured sampled flow rate and Api and Ap^ are 
measured pressure differences of the main flow stream 
with and without sampling operating. This assumes that 



as tine sole lonknown in the equation ^ 




where g 
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the density of the sampled flow qr is the same as the 
density of the total flow d. 



Thus Qi is calculated as t^i ^ 9 r— 71 h • From 

V^p./Apj -1 

this equation, the error 5Qi in Qx can be determined as 
« . However , as .^Apj/Ap^ — ^ 1 / - Therefore , in 

Qi A Qi 

order to avoid excessive errors in the measurement of the 
total flow rate Qi when using the method of WOOO/49370 it 
is important that the difference between Api and Ap2 is 

10 large, and consequently that the flow rate of the sampled 
flow is a substantial proportion (i.e. more than 5%) of 
the total flow rate. When such sampling proportions are 
used in high flow situations, for example in a high-rate 
gas condensate wells, they can overwhelih the separator, 

15 or alternatively require a larger, and therefore, more 
expensive separator. 

The present invention is at least partly based* on 
the realisation that isokinetic san^ling from a fluid 
stream can be achieved by monitoring the ratio of the 

20 flow x-ate of the san5>led flow to the total flow rate of 
the main fluid stream and adjusting the sampling rate 
accorclingly . Contrast this with WOOO/49370 which 
providies no suggestion that any repeated or continuous 
adjustment • to achieve isokinetic sampling should be 

25 carried out, and indeed assume that san^ling is carried 
out under isokinetic conditions in order to measure the 
total flow. rate. 
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Thus, a first aspect of the present invention 
provides a method of sairpling a multi-phase fluid stream 
comprising the steps of: 

saic^ling, with a sampling probe, a portion of the 
5 flu±d stream; 

measuring the flow rate Qs of said sampled portion; 

and 

measuring, independently of the san^ling step, the 
total flow rate Q of the fluid stream, 
10 wherein the flow rate of the sampled portion is 

controlled according to the ratio of the flow rate of the 
sampled portion to the flow rate of the fluid stream, in 
order to obtain substantially isokinetic san^ling of the 
fluid stream. 

15 Preferably the flow rate of the sanqpled portion is 

controlled such that the equation — = is satisfied 

to within 10%, more preferably to within 5%, where J^rote 
and Apipe are the respective cross-sectional areas of the 
sampling probe opening and the fluid stream at the 

20 sampling point . 

The flow rates Q and Qs may be mass flow rates, and 
the flow rate Qs of the sampled portion of the fluid 
stream is preferably less than 0.2% of the total flow 
rate Q, more preferably less than 0.05% of the total flow 

25 rate 0. 

The method may be used in high-rate gas condensate 
wells, where the gas flow rate is at least 5 MMscf/day. 
In wells with higher flow rates, the flow rate of the 
sampled portion of the fluid stream is particularly 
30 preferred to be a small proportion of the total flow rate 
so that downstream analytical processes on the sampled 
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portion are not overwhelmed. An advantage of tlie present 
inventi.on is that it can be used to san^le relatively 
small p>roportions of the total flow. 

Pr^eferably the method further includes the step of 
5 separating said sampled portion into a sampled liquid 
flow aixd a sampled gas flow, and said step of measuring 
the sampled flow rate includes measuring the flow rate of 
the sampled liquid flow and the flow rate of the sampled 
gas flow. Either or both of these sampling flows may be 

10 measured using a Coriolis flow meter. This allows more 
accurate measurement of the flow rates of the different 
phases of the fluid stream. A Coriolis flow meter also 
allows the density of the flow to be measured. 

In. certain embodiments the step of measuring the 

15 total flow rate includes measuring the density of the 

fluid stream. The step of measuring the total flow rate 
may include measuring a differential pressure, which may, 
for excLxnple, be measured using a Venturi meter, or be an 
impact pressure measured using a Pi tot tube. 

20 Pr^ef er ably . the step of measuring t^ie total flow rate 

is carried out continuously or repeatedly during the 
sanpling process . 

Preferably the step of controlling the flow rate of 
the sairepled portion is carried out continuously or 

25 repeatedly during the sampling process. 

The continuous or repetitive measuring of the total 
flow rate and controlling of the flow rate of the san^led 
portion allows correction of the flow rate of san^led 
portion in response to changes in the mass flow rate . of 

30 the main flow stream. 

In this respect, the term ~ continuous" also includes 
repeated discrete applications of the measuring and 
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controlling steps above at regular intervals, these 
intervals being shorter^ than the time scale over which 
significant changes in the total mass flow rate of the 
main stream occur* 

The above method is particularly applicable where 
the fluid stream is a feed from a well-head, and may be 
carried out at a well-head site. 

Preferably the method further includes, prior to 
sampling the fluid stream, the step of conditioning the 
fluid stream. 

Another aspect of the present invention provides a 
sampling apparatus for a multiphase fluid stream which 
can carry out the method of the above aspect. 

According to this aspect, the present invention 
provides a sampling apparatus for sampling a multi-phase 
fluid stream, con^rising: 

means for measuring the total flow rate Q of the 
fluid stream; 

a. sampling probe for sampling a portion of said 
fluid stream; 

ei meter for measuring the flow rate Qs of said 
sairpled portion; and 

a controller adapted to control the flow rate of the 
sampleij portion according to the ratio of the flow rate 
of the sampled portion to the total flow rate of the 
fluid stream, in order to obtain substantially isokinetic 
sampling of the fluid stream. 

Preferably the controller is adapted to control the 
flow rate of the sampled portion such that the equation 




is satisfied to within 10%, more preferably to 
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within 5%, where Aprobe and ^pipe a^re the respective cross- 
sectional, areas of the sampling probe opening and the 
fluid strream at the sampling point. 

Flow rates Q and Qs may be mass flow rates ^ in which 
case the meter is a mass flow rate meter. 

Preferably the cross-sectional area of the sampling 
probe opening is less than 0.2% of the cross-sectional 
area of the fluid stream at the sampling point, and thus 
the flow rate of the sainpled portion is a similar 
fraction of the total flow rate of the fluid stream. 
More preferably the cross-sectional area of the sanpling 
probe opening is less than 0.05% of the cross-sectional 
area of the fluid stream at the sampling point. 

The sampling apparatus may form part of a high-rate 
gas condensate well wherein the gas flow rate is at least 
5 MMscf/<iay. 

Preferably the apparatus further includes a 
separator- for separating the sampled portion into a 
sampled liquid flow and a sampled gas flow; a meter for 
measuring- the flow rate of the sampled liquid flow; and a 
meter for measuring the flow rate of the sampled gas 
flow. Either or both of these meters may be Coriolis 
flow meters. 

In some embodiments the means for measuring the 
total flow rate may include a densitometer for measuring 
the density of the fluid stream. The means for measuring 
the total flow rate may additionally or alternatively 
include a differential pressiire sensor, such as that for 
use with a Venturi meter or a Pitot tube. 

Preferably the means for measuring the total flow 
rate operates continuously or repeatedly during the 
sanpling x>^ocess. 
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Preferably the controller controls the flow rate of 
the sampled portion continuously or repeatedly during the 
sair^ling process. 

Preferably the sampling probe is made from damage 
resistant material such as timgsten carbide or cobalt 
alloy, so that the cross sectional area of the* saitpling 
probe opening does not vary during the sair5>ling process. 

The sampling probe is preferably located in the 
centre of the pipe section. 

Preferably the sampling apparatus further includes a 
flow coaditioner located upstream of the sait^ling probe. 

The fluid stream may be a feed from a well-head . 

In order to achieve representative isokinetic 
sampling of a multi^-phase stream using san^jling at a 
single sampling point, it is preferable that, at the 
designated sairpling cross-section of the pipe, the cross- 
pipe GLR profile is laniform. 

Previously flow- focus sing mixers have been used to 
try to provide such a uniform GLR profile. Exan^les of 
such inixers are plug-valve mixers, for exait^jle as shown 
in GB 2301297 and US 5894080. A similar mixer to that in 
GB 23 01297 is shown schematically in Figure 2a. Where 
similar reference nxomerals are used in Figtire 2a, they 
refer to features already described in association with 
Figure 1. 

However, use of flow- focussing mixers, such as the 
mixer 1015 shown in Figure 2b, can result in too much 
liquid being entrained into the central pipe region as a 
result of the focussing action of the mixers. At low 
pressures there may be excessive liquid jetting, leading 
to very non-uniform cross-pipe GLR and a large GLR error. 
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Figiare 2b shows the GLR cross-pipe profile for air/water 
flows with an equivalent CGR of -60 bbl/MMscf gas (gas 
rate qg ^ 1400 m^/h or about 1-2 MMscf /day) , and shows 
that the measured GLR has an error of approximately -70% 
5 from the reference value at the pipe centre. 

It is anticipated that as the wellhead pressure 
increases, the GLR cross-pipe non-uniformity will improve 
as a. results of the reduction in the gas- liquid density 
contrast, and hence in the gas-liquid velocity slip. 
10 However, the adverse effect of liquid jetting will still 
remain^ even at high wellhead pressures. Similar results 
are expected to be seen with the plug- valve type mixers 
referred to above. 

15 Therefore a f\irther aspect of the present invention 

provides a flow conditioner for a wellhead sampling 
system, the flow conditioner being a pipe section 
including: a swirl inducing section; a flow straightener; 
and an orifice plate. 

20 To achieve accurate isokinetic sairqpling at wellhead 

at a single sampling point, it is preferable to have a 
unifoann cross-pipe GLR profile at the designated sampling 
cross -section. For annular or annular-mist flows 
prevailing, in a wellhead sampling system (i.e. not fully 

25 developed vertical upward or downward flow) , this can be 
achieved by removing liquid film from the pipe wall and 
by entraining it homogeneously in the gas core. It is 
also preferable that liquid drop re-deposition onto the 
pipe, wall should be negligible at the sampling cross- 

30 section. Drop size is preferably small thereby 

mindLmizing gas-drop velocity slip and the rate of drop 
deposition. These factors will lead to a more 
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representative sampling of gas and liquid flow rates, and 
hence a better measure of the GLR. 

The creation of the initial swirl, upstream of the 
orifice plate, preferably creates a uniformly 
-5 circumfeorentially distributed liquid film on the wall of 
the pipe section. The orifice plate sheds any such 
liquid film on the pipe wail back into the gas flow. 
After the initial formation of a liquid rich core 
downstream of the orifice plate, t-urbulence causes good 

10 mixing of the licjuid and gas phases. Straightening the 
flow of tbe gas core prior to the orifice plate reduces 
the amount of liquid that is redeposited onto the pipe 
wall after shedding by the orifice plate. Therefore a 
relatively 'uniform cross-pipe GLR can be achieved, making 

15 it less important where on the pipe diameter in the 
sampling cross section a sampling probe is placed. 

The swirl inducing section may include a tangential 
input to the pipe section. Alternatively the swirl 
inducing section may include a helical insert inside the 

20 pipe section. 

The flow straightener is preferably located in the 
centre of the pipe section and such that, there is an 
annular g-ap between the flow straightener and the wall of 
the pipe section, thereby straightening the gas core of 

25 the fluid, stream whilst not affecting any liquid film 
located on the wall of the pipe. 

The orifice plate may have a diameter of its open 
section w^hich is between 0.85 and 0.95 of the local 
diameter of the pipe section, preferably between 0.9 and 

30 0.95 of tbe local diameter of the pipe section (also 
referred to as the diameter ratio /?) . 
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The pipe section may narrowed at the orifice plate. 
The. narrowed pipe section gives a higher gas flow 
momentxim and higher gas Weber nxjmber- 

To in^rove CGR measurement accuracy of lean gas 
5 condensate wells (CGR = 0 to 50 bbl/MMscf ) , a dual- 
orifice mixer arrangement is preferably used for reducing 
liquid rate under-sampling at low liquid loadings. 

In one exctrrtple of such an arrangement^ the first 
(i.e. upstream) orifice plate has a diameter ratio such 
10 as i8i=0. 9, and the second orifice plate has a diameter 
ratio such as ^2=0* 95. The two orifice plates may be 
axially separated by ID to 3D, where D is the basic 
diameter of the pipe at the orifice plates. 

The first, upstream orifice plate (diameter ratio fii) 
15 is intended to shed as much as lic[uid to the pipe centre 
as pdssilDle (with a relatively small p -ratio) ; whilst the 
second orifice plate is used to shed any relatively thin 
liquid film that might be caused by the re-deposition of 
the liquid drops shed by the first orifice plate, after 
20 travelling a few pipe diameters . 

A sampling probe may be used to sample the fluid 
stream, and is preferably located downstream of the 
orifice plate (s). More preferably the sampling probe is 
located at a distance between ID and 3D downstream of the 
25 downstream orifice plate, where D is the diameter of the 
pipe sect;ion at the orifice plate (not the diameter of 
the orifice plate itself) . 

Preferably the sampling probe is located in the 
centre of the pipe section. 



30 
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Anotiaer aspect of the present invention provides a 
method forr conditioning the flow of a multi-phase stream 
comprisinQT the steps of successively: 

forini_ng a liquid film at the wall of a pipe section 
which carr-ies the fluid stream; 

strai^ghtening the flow of at least a core of the 
fluid stream; and 

shedcLing the liqiiid film into said core. 

Prefesrahly the liquid film is substantially 
uniformly circiimf erentially distributed on the pipe wall. 

The step of shedding the liquid film may use at 
least one orifice plate. The or each orifice plate may 
have an open diameter of 0.85 to 0.95 of the diameter of 
the pipe section in which it is placed, preferably 
between 0. 9 and 0.95 of the diameter of the pipe section. 
The pipe section may narrow at at least the orifice 
plate. 

The step of forming the liquid film may include 
causing the fluid stream to swirl about the axis of the 
pipe section. Such rotary, motion may be achieved by, for 
example, tangentially introducing the fluid stream to the 
pipe section, or by passing the fluid stream through a 
helical insert inside the pipe section. 

Prefearably the step of straightening the core of the 
fluid stream does not straighten a portion of the fluid 
stream adjacent to the walls of the pipe, and therefore 
does not a£fect the liquid film. 

The fluid stream may be sait^led using a sampling 
probe, the san^ling preferably taking place downstream of 
the orifice plate (s) . More preferably the sampling probe 
is located at a distance between ID and 3D downstream of 
the downstream orifice plate, where D is the diameter of 
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the pipe section at the orifice plate (not the diameter 
of the orxfice plate itself) . 

Preferably the sait^ling probe is located in the 
centre of the pipe section. 
5 The flow conditioner or method of the above two 

aspects may be used with the sampling methods and 
apparatuses of the earlier aspects to improve the 
acciiracy of the isokinetic sampling. 

10 • - Embocdiments of the present invention will now be 

described xn relation to the accompanying drawings, in 
which : 

Figure 1 is a schematic illustration-' of a known 
wellhead sampling system, and has already been described; 
15 Figure .2a is a schematic illustration of a flow- 

focusing mixer similar to that shown in GB 2301297 the 
prior art and has already been described; 

Figure 2b shows the GLR cross-pipe profiles 
downstream of the mixer of Figure 2a and downstream of a 
20. flow conditioner according to an embodiment of the 
present invention; 

Figure 3 is a flow chart showing a method according 
to an embodiment of the present invention; 

Figure 4 is a schematic illustration of an 
25 isokinetic satmpling system according to an embodiment of 
the present invention; 

Figure. 5 is a schematic illustration of an 
isokinetic sait^ling system according to another 
embodiment of the present invention; 
30 Figure 6 is a schematic illustration of an 

isokinetic sampling system according to a further 
emboddLment. of the present invention; 
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Figure 7 is a schematic diagram of a flow 
conditionear according to an embodiment of the invention; 

Figure 8 is a schematic diagram of a flow 
conditioner- according to an alternative embodiment of the 
invent ion ; and 

Figure 9 is a schematic diagram of a flow 
straightener- which may be used in the flow conditioners 
shown in Fignres 7 and 8 . 

Elements of a method of sampling according to an 
embodiment of the present invention are shown in the flow 
chart of Figure 3- 

A portion of a fluid stream is sampled, and the flow 
rate Qs of that sampled portion measured. The total flow 
rate Q of the fluid stream is also independently 
measured, and the ratio of the two measured flow rates 
calculated. Based on this ratio, the sampling flow rate 
(i.e. the proportion of the fluid stream which is 
sampled) is adjusted in order to substantially obtain 
isokinetic sampling conditions. This process is 
preferably either repeated or carried out continuously 
throughout the sampling process. 

The sampled portion may be analysed concurrently 
with the flow rate measurements and adjustments, for 
example to provide the GOR, the WLR and samples for P\n? 
analysis. This analysis may also use the flow rate 
measurement of the san^led portion. 

Figure 4 shows an isokinetic sampling system 
according to an embodiment of the present invention. 

Ah isokinetic sampling model is used to calculate 
the target sampling flow rate Qso. according to the main 
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stream, flow rate 0©^^ . One example of such an isokinetic 
sarnpling moidel is -^^= ^^* for substantially uniform 
mixed flow." 

In contrast to the method described in relation .to 
Figure 1 above, in the present embodiment, the main 
stream flow rate Q^tt is measured on-line in the pipe by a 
main pipe flow meter 40. This may be, for example, a 
Venturi or an orifice meter. The measured Q^ff is then 
fed into the sampling model so that Qso is updated as Q^^^ 
varies. 

When the on/off valve 35, allowing sampling to take 
place, the actxial sampling flow rate Qs is measured by a 
san^ling flow meter 30 positioned in the sampling line. 
The actual samipling flow rate Qs and the target sanpling 
flow rate Qso are compared in comparator 33, which creates 
a flow control signal 31, on the basis of which control 
valve 32 is adjusted by the actuator 34 so as to keep Qs = 

Alternatively flow control signal 31 may be created 
as a result of a direct comparison between the main 
stream flow rate Q^ft measured by the main pipe flow meter 
40 and the sampling flow rate Qg measured by the sampling 
flow meter. 

The main pipe flow meter 40 of the embodiment of 
Figure 4 may be positioned upstream or downstream of a 
mixer or other flow conditioner (not shown) . In certain 
eznbodiments, tlae main pipe flow meter will also serve as 
a good gas-liqiaid mixer, and so a separate flow 
conditioner may not be required. 

However, u.se of a flow conditioner upstream of a 
flow meter 40 can searve to eliminate or substantially 
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reduce the gas-liquid velocity slippage for the flow 
passing through the main pipe flow meter 40, thereby- 
allowing a Inomogenous flow model (e.g. steuadard single- 
phase Venturri equations) rather than a slip -model (e.g. a 
5 2 -layer, separated model for annular or annular-mist 

flows) to be used to calculate the main stream flow rate. 
The homogenous mixing state is assumed to be maintained 
over some distance downstream of the main pipe flow meter 
40, enabling the sampling probe to perform the desired 
10 isokinetic sampling. 

Prefer^red embodiments of separate mixers or flow 
conditioners are described below. 

Another embodiment of the present invention is shown 

15 scheraaticallLy in Figure 5 . An embodiment of the method 

of the present invention is also described in relation to 
this Figure . Again components which have already been 
discussed in relation to Figrures 1 and 4 have the same 
reference nuunbers . 

20 This embodiment has a flow meter 40 which combines a 

Venturi pressure drop measurement 41 with a gamma-ray 
densitometer 42. An example of such a meter is the* 
PhaseTester/PhaseWatcher (TM) Vx Venturi meter described 
in SPE paper 77405 (by S. Jayawardane and B.C. Theuveny 

25 ''PVT Sampling With Multiphase Flowmeters-Theoretical 
Justifications and Field Limitations", SPE Annual 
Technical Conference & Exhibition, San Antonio, Texas, 29 
September to 2 October 2002) . The pressure euid 
temperature (p, T) of main flow stream 11 can also be 

30 measured at this point. The Vx flow meter has been shown 
(see SPE paper 76766) to give better multiphase flow rate 
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measureznenti than a test separator, especially at high gas 
rate and/or unsteady flow conditions. 

This embodiment uses a total mass flow rate 
measurement Q = K^y[^p^p^ , obtained from the combined 
measurements of high-energy mixture density (p^, 
calculated . to about 10 kg/m^ accuracy by the Vx nuclear 
source) and Ap from the Venturi measurement (calculated 
to within 1% accuracy) to set the extraction rate for 
isokinetdLc sampling, where K is a calibration factor 
depending on Venturi throat diameter, inlet diameter and 
the discfciarge coefficient; see, for example, eq.5 in SPE 

paper 774:05. The mass flow rate error ^ = if A(Ap) ^ ^Pg] 

is thus within 2.5% foir a mixttire density pn = 300 kg/in^. 

The sampling probe 2 0 has an opening area (Aprobe) . or 
ratio of its area to the pipe area (Apipe) , which is known 
accuratel-y to within 2% and is maintained unchanged over 
a long period- To achieve this, the probe tip 18 is made 
from dama.ge-resistant material such as tiingstea carbide 
or cobalt, alloy. 

Consideration is also given to the size of the 
opening of the sampling probe tip 18. In order to keep 
the relative error in the cirea of the opening of the 
probe, air in the area ratio to within 2%, a 1% accuracy 
is requir^ed in the diameter of the opening of san^ling 
probe 20 . Therefore for a small opening diameter 
sanpling probe 20 (e.g. 1.5 mm), the tolerance in 
diameter lias to be very high {< 15 um) . Increasing the 
opening d.iameter of the sampling probe 20, for example to 
3 mm, red-uces the demand on the machining tolerance for 
the probe (to < 30 |im for the 3 mm example above) . For 
high GLR flows, a larger probe size may also be 
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beneficial as it will increase the amount of liquid 
sangpled and therefore ease the sizing of the Coriolis 
mass flow meters 30a and 30b for use at higher ranges of 
gas anca licguid flow rates in the mini-separator 50. 

It is preferable to have a uniformly mixed wet-gas 
flow at: the sainpling cross section. This is achieved by 
placing a flow conditioner 2 00 upstream of sampling probe 
20. Pr^ef erred embodiments of the flow conditioner 200 
are described below. 

It is also preferable that the stream at the Venturi 
portion of the Vx flow meter 4 0 is also homogenous wet- 
gas flow (i.e. with negligible gas-liquid slip), and so a 
further^ flow conditioner may be provided upstream of the 
Vx flow meter 40. 

With the on/ off valve 3 5 open and operating in 
conditions of well-mixed wet-gas flows in which there is 
low Xicgoiid mass (i.e. Ql/Qg < 1) the total mass flow rate 
measurement Q from the Vx flow meter 40 

{Q = K^^.p„ ^Qa'^QL^Qci^ + QLlQc) ) is first used to set an 
approximate ^dry' gas sample extraction masis flow rate 

by adjusting flow control valve 32 in 
the mini -separator 50 through use of the actuator 34. 

The mass flow rate of the 'dry' gas sample flow Q^^, 
after gas-droplet separation in the gas-drop separator 
51, can be directly measured using a Coriolis flow meter 
30a. 

The sampled liquid mass flow rate Q^j^ is measured 
using an instrumented liquid collection system 52 in the 
mini -separator 50. Precise achievement of isokinetic 
samplingr has less effect on the sample liquid mass flow 
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rate Qs^l than on the sample gas flow rate because of the 
higher inertia of liquid droplets compared to the gas . 

In the mini^separator 50 -the sampled liquid mass 
flow rat:e Qs,l can be derived from the change in the 
liquid level oyer time b(t) in the liquid collection 
vessel 52 from a Dp level sensor (not shown) 

) r where Aj, ±s the cross-sectional 
area of the liquid collection vessel 52) . 

However, in the embodiment illustrated din Figure 5, 
is also measured using a second Coriolis mass flow 
meter 3 Ob, which also gives the liquid mixture density 
Ps,L' For example the second Coriolis mass flow meter 30b 
may be sl Promass 80A« Coriolis meter as made by Endress- 
Hauser. This meter can measure liquid mass flowrate to 
±0.1%, gras mass flowrate to ±0.5% and density to ±1 kg/m^, 
and is reported to measure flow rates as small as 0.1 
kg/h. 

An extra linear array of MEMS (Micro Electro 
Mechanical Systems) densitometer can be installed along 
the heigrht of the liquid collection vessel 52, allowing 
oil and water densities ips,o, Ps^w) to be determined from 
the sepa.rated oil-water interface profile at separator 
,temperat:ure and presstxre (Ts, ps) . 

It is preferable that the liquid is separated at 
line pressure not following a pressure reduction to 
minimize mass transfer between the phases. It is also 
preferable that the lic[uid sanqple mass flow rate 
measurement system, whether it is a Coriolis meter or a 
measurement bottle has an error of 2% or less to achieve 
good results in the closed- loop feedback system described 
below (a. better than 2% accuracy is necessary to obtain 
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an accurate measurement of the sampled liquid flow rate 
Qs,L* and of the main-pipe liquid rate Gl =(a^ /A^^jx^^^ ) . 

The wetter cut (yvf) can then be determined from the . 
liquid mixture density relation =?^wA,w +(l~yw)/^.s.o - 
5 Alternatively the water cut can be measured by other 
means, for example using microwave or optical sensors, 
such as the RedEye™ infrared water- cut meter from 
eProduction Solutions, Texas US, prior to complete 
separation of the liq[uids. 

10 The f ±rst measure of the liquid sample mass flow 

rate j25.t is then combined with the measured sample gas 
mass flow rate in the controller 60 to obtain the total 
sample mass flow rate - The controller 60 then 
compares this with the current main pipe total mass flow 

15 rate as measured by the main pipe flow meter 40. The 

sample ext^raction rate 19 is then varied as necessary by 
adjusting the flow control valve 32 in the mini- separator 
50, througli use of the actuator 34, so as to cause the 
sample mass flow rate to approach {Ap^/A^)xQ, and thus 

20 isokinetic sampling conditions. 

The above process can then be repeated for 
successive measurements of the total sample mass flow 
rate, thereby achieving convergence with the isokinetic 
sampling conditions. The process may be run effectively 

25 continuously throughout the sampling process, and thereby 
adjust the sampling flow rate as necessary when the flow 
of the main stream varies for any reason. In this 
respect, the term continuous" includes repeated discrete 
applications of the measuring and correcting steps above 

30 at regular intervals, these intervals being shorter than 
the time scale over which significant changes in the 
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total mass- flow rate of the main stream occur. A typical 
interval is less than 30 sec, and preferably less than 10 
or 5 sec. 

For the maintenance of isokinetic san^jling in this 
system, the existence of a pressure difference between 
the main pipe 220 and the mini -separator 50 does not 
affect the sampling. Mass transfer between the liquid 
and gas phases can occur without prejudicing the sampling 
control process since the process is based on total mass 
flows and not on the volume rate of a single phase. 

However, it is preferable that the pressures in the 
mini -separator 50 and in the main pipe 220 are 
approximately equal, so that no mass transfers between 
the gas and liquid phases occur between the pipe and the 
separalior. ■ In this situation, the main pipe gas and 
liquid mass flow rates can be obtained as follows: 

Total mass flow rate Q^K^lSp.p^ , obtained from the 
Vx Venturi Ap measurement and the high-energy mixture 
density measurement; 

-L.±quidmass flow rate QL^(\ipef ^rote^QsM ^ obtained 
from the isokinetic sampling; 

Gas mass flow rate Qo-^Q-QL^Q-{A,ipJ A^)^Qsm' 
preferably obtained from the Vx total flow measurement 
and isokinetic sairpling, also obtainable purely from the 
isokinetic sampling as Qc^[\ipef \rate)xQsj^ > 

which can 

therefore be used as a quality control check. 

As a result, the PVT relations are no longer 
required to convert the liquid mass flow rate measured at 
separator conditions (ps/ Ts) to that at main stream 
conditions (p, T) in order to obtain gas and liquid mass 
flow raties from the Vx measurements and interpretation 
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models. However, PVT relations are reqpzired to convert 
the mass flow rates at line conditions (p, T) to volvune 
flow rates at standard conditions. 

The mini -separator 50 also provides the following 
5 measurements of the sanpled flow: 

Densities Ps,g/ Ps,l, Ps.o. Ps.w (of gas, iiqijid, oil 
and water samples respectively) obtained from the 
Coriolis gas and liquid flow meters 30a and 30b and the 
MEMS densitometer array reading from the licjuid 
10 collection vessel 52; 

Water cut Vw ={Ps,l- Ps,oV{ps,w - Ps.o) calculated from the 
densities above; 

GLR = {Qs,g^ Ps,g)KQs,l^ Ps^) calculated from the separator 
gas and liquid vol\ime flow rates; 
15 Iiiquid mass loading m-Qsj^/Qs^ calculated from the 

separator flow rates - ideally m < 1 for a well mixed 
wet— gas regime which the system is being used to analyse ; 
and 

PVT data from the scunples obtained from the 
20 separated liquid flow and separated (i.e. ^dry') gas flow 
in the mini -separator 50. 

Another embodiment of the apparatus and method of 
the present invention is illustrated schematically in and 
25 witli reference to Figure 6. 

As with the embodiments described above, it is 
preferable that the wet-gas flow at the sampling cross- 
section is well-mixed, and this is achieved using mixer 
or flow conditioner 200 positioned upstream of the 
30 samfE^ling cross-section. Again, preferred embodiments of 
the flow conditioner are described below. 
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An averaging Pi tot -tube (APT) 44 is used to derive 
the main pipe total mass flow rate when operating at 
xmpact -pressure measuring (Pitot) mode, as well as as a 
san^^ling probe. The APT 44 has multiple holes, for 
example four or six holes, specifically arranged along 
the diametrically aligned tubing. For exaic^jle, for a 
six-hole probe the holes are preferably placed on the 
probe txibing at the following positions: -.913R, -.707R, 
-.408R, +.408R, +.707R, +.913R, where R is the radius of 
the pipe. The hole positions preferably correspond to 
the centres of the sait^led equal area annuli of the pipe 
circular cross section. As with the previous embodiment, 
for sampling accuracy it is preferable that the area of 
the hole openings of the APT 44 (Aprobe = Ntoie x A^oie) , 
relative to the pipe area Apipe at the sampling cross- 
section, is known to within 2% and does not change over 
the time of testing and sampling. This can be achieved 
by making the APT 44 from damage resistant material such 
as tungsten carbide or cobalt alloy. 

One exan5>le of an APT is the Rosemount Probar^ 
insertion meter operating at impact-pressxire measuring 
(Pitot) mode only, which claims a relative accuracy of 
±1.1% in voliametric flow rate measurements for steam, 
liquid and gas. APTs can also be obtained from 
Thermocouple Instruments Ltd. Of Cardiff, UK. 

An APT, probe has very low pressure loss and 
generally costs less than a Venturi flow meter. The pipe 
averaging characteristics of the APT 44 reduce any 
variations in the cross-pipe flow at the sampling cross 
section. 

A high energy gamma-ray densitometer 43 is installed 
upstream of the APT probe 44 and downstream of the flow 
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conditioner 200 to obtain the density p^. The 
densitometer 43 can have an accuracy within 10 kg/m^ from 
a 30 mCi ^^^Cs source (with high energy gamma-ray at 
661 .66 KeV) and a Nal scintillation detector . 
5 The APT 44 is first set to operate in normal impact 

P2ressure measurement mode, with on/off valve 35 set in 
ttie off position and with no saitpling flow being taken. 
Thie APT measures the differential pressure Ap across the 
main-pipe static pressure tap and the APT probe static 
10 prressure tap (typically to within 1% accuracy) , and the < 
main pipe total mass flow rate is calculated as 

Q = K^^Ap.Pff (here Ki is a calibration factor dependent on 
thie probe geometry and flow Reynolds -number) . 

The APT 44 is then switched to operate in sampling 

15 mode by opening the on/off valve 35 and by adjusting flow 
control valve 32 through actuator 34, The averaging 
nature of the APT 44 means that the main flow stream is 
sampled at multiple points across the diameter of the 
san^ling cross section. Therefore variations in the 

20 crross-pipe GLR profile, which may be fo\xnd downstream of 
an imperfect flow conditioner, can be largely overcome. 
Alternatively, when an APT is used, as in the current 
embodiment, the performance requirements of the flow 
conditioner can be made less stringent. 

25 Operating in conditions of wet-gas flows in which 

thiere is low liquid mass (i.e. Qc/Qg < 1) the total mass 
flow rate measurement Q from the gamma-ray densitometer 
43 and the APT 44 operating in impact pressure 

measurement mode ( K^^Ap.p„ -Qg-^Ql-Qg^^QlIQg) ) is 
30 fi-rst used to set an approximate 'dry' gas sample 
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extiraction mass flow rate Qsjs^ "^{.^rvbe^ \qK)^Q adjusting 
flow control valve 32 in the mini -separator 50. 

Tiie remainder of the san^ling process is the same as 
that described in relation to Figure 5 above. In 
5 particTilcir the flow control valve 32 is adjusted via the 
actuator 34 so that the sanqpling flow 19 is substantially 
isokinetic with the main pipe flow stream 11. 

In a further embodiment, a single-hole sampling 
probe 20, such as that used in the embodiment described 

10 in relation to Figure 5 above, can also be set to operate 
in ixnpact pressure mode by switching valve 35 off and 
measuring the differential pressure hp across a main pipe 
static pressure tap and a static pressxire tap on the 
probe 20 , This measurement of 23p could then be used to 

15. derive an independent estimate of the main pipe total 

flow rate using the density Pb measured by. densitometer 
42. 

In a further embodiment, in applications where the 
probe pressure loss (i.e. from the probe opening to the 

20 probe pressure tap) is relatively small, the isokinetic 
saxigplitig condition can be substantially achieved by 
nullxng the hp described above by adjusting the sarqpled 
flow rate of the saonpled portion through the use of flow 
Gontzrol valve 32 « In this embodiment, the Ap signal can 

25 be used as a control signal to control acttiator 34, thus 
avoic3ing the use of f low meters 30a and 30b- Note that 
this exiibodiinent can be adapted to either the single hole 
probe shown in Figure 5 or the zmaltiple hole pi tot tube 
probe shown in Figure 6. 

30 Note that the multiple-hole pitot tube sampling 

embodiments described in connection with Figiire 6 can 
also be used in applications where the independent 



wo 2005/031311 



PCT/GB2004/003849 



29 . 

measurement is made using Vx. type meters as described in 
connection with Figrure 5- 

Figure 7 is a diagram of a pipe section which is 
aroranged to provide a relatively \iniform cross-pipe 
profile at the sampling cross-section, particularly with 
regard to the GLR. The pipe section shown has a flow 
conditioner 2 00 for providing such a cross-section. 

The flow conditioner 200 has three elements: a swirl 
incducing portion, which is provided in this embodiment by 
the tangential flow input 201; a gas flow core 
straightener 202; and an orifice plate 203. 

The tangential flow input 2 01 causes the entire flow 
to swirl, thereby creating a substantially uniformly 
circumf erentially distributed liquid film on the pipe 
wall 205. This removes variations due. to the inflow 
conditions, such as non-axi symmetric gas and liquid flow 
di s fcr ibut i ons'. 

The core of the flow, which is mostly gas, is then 
straightened by the core straightener 2 02 positioned 
downstream of the tangential input 201- An example of 
such a core straightener 202 is illustrated in Figure 9 
anca. comprises a conventional straightener 212 (such as 
those described in US 3840051, US 5309946 and US 5762107) 
positioned with an annular gap 213 between the 
straightener 212 and the pipe wall 205 so as not to 
interfere with the liquid flow along the pipe wall- The 
core straightener 202 removes the swirl in the central 
gas core of the flow, thereby reducing the amount of 
•licjuid re-deposited on the wall after the orifice plate 
mixing- This facilitates the obtaining of consistent 
gas- liquid flux profiles at the sampling cross-section. 
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Downstream of the core straightener 202, an orifice 
plate 20 3 is positioned, which causes the liquid film 
flowing on the pipe wall to be shed back into the now 
straightened gas core. The dicimeter ^-ratio (the ratio 
5 of the open diameter of the orifice plate 203 to the 
diameter of the pipe) is preferably between 0.85 and 
0-95. A. further orifice plate (not shown) may be 
positioned downstream of orifice plate 203, and sheds any 
liquid film re-deposited on the wall after the orifice 

10 plate 20 3. 

The liqu.id shed from the pipe wall by the orifice 
plate 20 3 is then subjected to an acceleration-induced 
catastrophic break-up by the gas flow which has a high 
axial velocity. In the region immediately downstream of 

15 the orifice plate 203,. a liquid rich cone is formed at 
the centre of the flow. Downstream of this cone, the 
turbulen^t kinetic energy is high and hence the secondary . 
turbulent break-up controls the droplet size 
distribiition. Further downstream, as the turbulent 

20 kinetic energy decays quickly, turbulent break-up ends 
and coalescence effects determine the liciuid droplet 
size. T*hus the cross-pipe liquid flxix distribution 
downstream of the orifice plate .203 is totally different 
from tha.t upstream of the plate 203. 

25 The sampling probe is preferably placed at an axial 

distance L downstream of the orifice plate 203, the 
distance X» being chosen so as to avoid sampling at the 
liquid-r'ich core region caused by the shedding of the 
liquid film t>y the orifice plate 203, and to allow for 

30 some break-up due to turbulence and. diffusion, leading to 
smaller drop sizes and more uniform cross-pipe GLR. The 
. preferred values of L are ID to 3D. 
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The use of an orifice-plate type mixer 203 is 
preferable to the flow- focusing type mixers, such as 
those disclosed in GB 2301297 and US 5894080. An easy- 
to— machine orifice-plate mixer 203 serves to shed the 
5 licjuid film from the pipe wall, which is then subjected 
to acceleration- induced catastrophic break-up by the 
incoming high axial- velocity gas flow. For a 4" pipe 
diameter (D = 100 mm) down-flow exaimple at low pressures, 
a much more uniform cross -pipe GLR has been 

10 experimentally measured 3D downstream of a single 

o*r±f ice-plate 203 with diameter ratio p= 0.95 (see Figure 
2b) than' that obtained with the flow- focusing mixer 15 
shown in Figure 2a (which is similar to that described in 
GB 2301297).. As a result, the GLR error at the pipe 

15 centre is reduced by about five- fold over the flow- 
focusing type mixer shown in Figure 2a. The orifice 
plate and the flow-focusing mixer were both tested at the 
same CGR of 60 bbl/MMscf, at a gas rate of 1.2 MMscf/day. 
The embodiment of Figure 7 shows the flow 

20 conditioner and sampling probe 20 implCTiented in a pipe 
section 210 of smaller diameter than the main pipe 220 
(for example a section of 2 inch (D = 50 ram) rather than 
4 Inch (100 mm) diameter) . For a given mass flow rate of 
a well stream, deploying the sanpling probe 20 at a small 

25 diameter pipe section 210 achieves the highest gas flow 
momentum (pcDb^) and hence the highest gas Weber number 
WBjd = PgUg^D/o (where o is the surface tension) . This may 
be beneficial as it: increases the shear for the liquid- 
film primary break-up from the wall, hence increasing 

30 drop formation; causes generation of small size droplets, 
thereby reducing the drop Stokes ntimber audi gas-drop 
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velocity slip; increases the cross-pipe unifonaity and 
the percentage of liquid entrainment . 

An alternative embodiment of the flow conditioner 
200 is shovm in Figure 8. This embodiment .is the same as 
that in Figure 6, except that the tangential input 201 
has been replaced by a bend inlet 221 followed by a 
helical insert 211 to provide the initial swirl. 

Alternative embodiments which cause axisymmetric 
distaribution of the lic[uid film on the pipe wall, by 
causing swirl in other ways may also be used. 

The above embodiments may form part of an overall 
sairpiing apparatus, being combined with the apparatuses 
and methods of earlier embodiments to provide a uniform 
cross-pipe GLR for sampling. 

While the invention has been described in 
conj-unctibn with the exerr^lary embodiments described 
above, many equivalent modifications and variations will 
be a:pparent to those skilled in the art when given this 
discILosure. Accordingly, the exenplary embodiments of 
the invention set forth above are considered to be 
illustrative and not limiting. Various changes to the 
descaribed embodiments may be made without departing from 
the spirit and scope of the invention. 



